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ABSTRACT 

We have discovered a class of eccentric binary systems within the Kepler data archive that have 
dynamic tidal distortions and tidally-induced pulsations. Each has a uniquely shaped light curve 
that is characterized by periodic brightening or variability at time scales of 4-20 days, frequently 
accompanied by shorter period oscillations. We can explain the dominant features of the entire class 
with orbitally-varying tidal forces that occur in close, eccentric binary systems. The large variety 
of light curve shapes arises from viewing systems at different angles. This hypothesis is supported 
by spectroscopic radial velocity measurements for five systems, each showing evidence of being in 
an eccentric binary system. Prior to the discovery of these 17 new systems, only four stars, where 
KOI-54 is the best example, were known to have evidence of these dynamic tides and tidally-induced 
oscillations. We perform preliminary fits to the light curves and radial velocity data, present the 
overall properties of this class and discuss the work required to accurately model these systems. 



1. INTRODUCTION 

Kepler, in its search for extrasolar planets, obtains 
high-precisio n time series phot ometry of 115.6 square de- 
grees of sky (jKoch et al.ll201C)l ) . With this new combina- 
tion of photometric quality and uninterrupted coverage 
of more than 160 000 stars, Kepler is discovering variabil- 
ity never before identified in stars. Here we present light 
curves that demonstrate a periodic brightening or vari- 
ability on time scales of days, frequently accompanied by 
harmonic oscillations at shorter periods. 

These light curves are remini scent of the t i dally dis- 
torted light curves theorized by iKumar et al.l (|1995[ ) for 
eccentric orbits. KOL54 (KIC 8112007), a highly ec- 
centric binary system composed of two A-type stars ob- 
served by Kepler (Welsh et al. 2011), was the first system 
discovered whose light curve is predominantly character- 
ized by these dynamic distortions. Binary tidal forces 
also induce p ulsations at intege r harmonics of the or- 
bital period (jKumar et al.l 119951: iWillemsl |2003[) : mak- 
ing KOI-5 4 the best example o f pulsations dr i ven in this 
way. See [H andler et al.' ()2002D . lDe Cat et al.l (2000) and 
iMaceromet al.. (2009. ') for other examples of systems with 
evidence of this type of driving. The discovery of KOI- 
54 has provided critical observational data to constrain 
the theories of dynamic tidal forces on stellar binaries 
(jFuUer fc Lail 120111: iBurkart et aIll2012D . but more ex- 
amples of these stars are required to generally test these 
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theories. 

Despite having unique light curves, we demonstrate 
that the 17 systems described in this paper are all tidally 
interacting binary systems in eccentric orbits. In 121 we 
describe the Kepler photometric time series of these sys- 
tems and propose a simple model to describe their dom- 
inant features in ij3l We fit our simple model to these 
systems and estimate their orbital param eters in ij3.1l 
demonstrating a technique proposed by IKumar et al.l 
(1995). In 2] we present spectroscopy for each system 
including radial velocity measurements of five systems 
that confirm they are in eccentric binary systems. Fi- 
nally we discuss the properties of this class of stars, the 
effects of the tidal interactions, and the remaining work 
to accurately model these systems in ^ 

2. KEPLER LIGHT CURVES 

We have discovered 17 objects in the Kepler data 
archive that demonstrate unique periodic variability. We 
refer to these objects as "heartbeat" stars because the 
shape of the variability for several are reminiscent of an 
echocardiogram. We show example light curves in Fig- 
ure [T] that exemplify the sudden and periodic variation 
in brightness. Notice that the heartbeat profile comes in 
many varieties and appears to be correlated with other 
periodic oscillations in the light curve of the star. 

For each system we analyzed a subset of the long ca- 
dence (one exposure taken every 29.4 min) light curves 
observed during Kepler quarters through 7 (MJD span- 
ning 54953.03 to 55 552.06). Each quarter is one fourth 
of the Kepler heliocentric orbital period of 372 d, except 
for the quarters and 1 (with lengths of approximately 
10 d and 33d respectively). Each subset contains enough 
data to constrain the heartbeat period to better than two 
minutes. See Table [T] for a list of long cadence quarters 
used to study each object. 

We used the light curves create d by the Kepler sim- 
ple aperture photometry module ("Twic ken et al.l 120101 : 
iR-aauclli fc Thomp son, 20 11.). These light curves are the 
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Figure 1. Example light curves of heartbeat stars found in the Kepler data. 



result of performing simple aperture phot ometry on the 
calibr ated pixels in the optimal aperture (jBrvson et alj 
120101 ). For more details on the calibration steps per- 
formed by the Kepler pipel ine see I Jenkins et a l. (2010), 
ICaldwell et all (poTof ). and lFanelhet al.l ()2011[ ). To re- 
move the long period effects that result from differen- 
tial velocit y abberation and the rmal effects from the 
spacecraft |Van Cleve et al.|[2011|) . we removed a piece- 
wise po lynomial fit from each quarte r using the program 
WQED (Thom pson fc Mullallvl 120091 ). We did not ad- 
just the lightcurves for contamination from other stars. 
The contamination is small (typically < 5 per cent) 
for these stars, as reporte d by the Kepler archive 
(|Fraauelli fc Thompsonll2011[ ). The contamination will 
have a small effect on the amplitude of the features we 
see in the light curves, but will not affect the shapes. 

We measured the heartbeat period by performing a 
nonlinear least-squares fit to the "comb" of equally 
spaced frequencies found in the power spectrum of each 
light curve. This comb results from the periodic, yet sud- 
den, changes in the brightness and from any harmonic 
pulsations present on the star. We fitted a sine wave at 
the orbital frequency found in the power spectrum as well 
as at a minimum of six harmonics of the heartbeat fre- 
quency to e stablish the heartbeat period (using Period04, 
iLenz fc Br eger 2005) (See Figure [S]). In performing the 
fit, the amplitude and phase of each harmonic were al- 
lowed to vary, while the frequency was fixed to an integer 
multiple of the orbital frequency. The fitted orbital pe- 
riods are reported in Table [1] At this stage we do not 
attempt to determine what fraction of the fitted ampli- 
tudes is due to the heartbeat shape and what fraction 
is due to stellar pulsation; it is not required to establish 
the period of the heartbeat. 

We folded each light curve at the measured period and 
then smoothed by performing a running average span- 



ning 15 minutes (see Figures |3] and U) . We set the zero 
phase to either the maximum or the minimum of the 
hearbeat feature, depending on the shape of the fea- 
ture. Because of the quality and length of the Kepler 
data, almost all the features shown in the folded light 
curves are real. However, because these are folded and 
smoothed curves, any variations that are not harmonics 
of the heartbeat period are lost. 

Notice that each folded light curve has a unique shape, 
quite different from KOI-54. Some dim before they 
brighten, others dim after they brighten, and others show 
distinct 'W' or 'M' shapes. Besides the heartbeat, stel- 
lar oscillations at harmonics of the heartbeat period are 
also apparent. For examp le, KIC 3230227, identified by 
lUvtterhoeven et al.l ()2011D as an eclipsing (5 Set star, has 
short period oscillations evident in the folded light curve 
and thus harmonics of the orbital frequency. Similarly, 
Figures [3] and S] reveal that most of the sample have pul- 
sations at harmonics of the heartbeat frequency. 

The presence of harmonic oscillations i s a clue to 
the na ture of these systems. For KOI-54, iWelsh et al.l 
(j2011| ) show that short period modes that are harmon- 
ics of the orbital frequency dominate the light curve 
once the effects of the tidal distortion and external heat- 
ing have been removed. Another hint to the nature of 
these objects comes from KIC 3230227, KIC 4936180, 
and KIC 8719324. They all show evidence of a graz- 
ing eclipse. Therefore, it is likely that the heartbeat 
shape is also caused by a binary companion. While sev- 
eral of the heartbeat systems identified h ere are also in 
the Kepler eclipsing b inary star catalog (jSlawson et al.l 
120111: iPrsaet all 1201 1[ ). all of those in the catalog (ex- 
cept KIC 3230227) show no evidence of an eclipse. 

In the next section we explain how this varied collec- 
tion of light curve shapes can be created with tidal dis- 
tortions in an eccentric binary system. 
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Figure 2. The amplitude spectrum of tlie Fourier transform of KIC 5034333. Eight evenly spaced harmonics are labeled as red lines at 
the top of the figure. To determine the heartbeat period, we fit the light curve at the fundamental frequency and this series of harmonics. 



3. TIDAL DISTORTIONS IN AN ECCENTRIC ORBIT 

Tidal forces ch ange the shap e of a star causing j p right- 
ness variatio ns (jBrown et all 1201 It iWelsh et all 120101: 
lMorrisiri985f) . For circular orbits, this results in a non- 
spherical star producing brightness variations at twice 
the orbital period, known as ellipsoidal variations. For 
eccentric binaries, the tidal force is phase dependent and 
the largest variations in shape and brightness occur at 
periastron. As a result, instead of seeing the double- wave 
in the flux variations, as in the case of the circular or- 
bit, the flux changes drastically when the binary system 
nears periastron. 

As a star passes through periastron, in the static tide 
limit, the star's shape adjusts to the instantaneous tidal 
force. The shape of the star's distortion can be ap- 
proximated by summing over axisymmetric and sectoral 
spherical degree two (^=2) modes. When observed at 
different inclinations, the relative contributions of the 
^=2, 171= —2, 0, 2 modes, integrated over the observable 
sphere, will change. Additionally, the observed phase of 
these modes will depend on the angle of periastron for 
the binary system. As a result the shape of the flux vari- 
ations due to the tidal distortions will depe nd on inclina- 
tion, angle of periastron, and eccentricity (jKumar et al.l 
[1995). The amplitude of these variations will depend on 
factors that influence the strength of the tidal forces, such 
as the masses of the objects and the orbital separation 
at peristron. 

The tidal brigh tness variations were derived by 
iKumar et al.l (|1995( ) while predicting the tidal oscilla- 
tions induced in a B-type star orbiting a neutron star. 
Even though the effect was too subtle to observe at the 
time, they did describe an analytic model whose shape 
depends on the orbital inclination, eccentricity and an- 
gle of periastron of the binary orbit. We calculate a grid 
of solutions to this model for different inclinations and 
periastron angles and plot each solution in Figure [5] for 
two different eccentricities. 

A side-by-side comparison of the grid in Figure [S] to 
the folded light curves in Figures |3] and |4] shows that this 
fairly simple model can reproduce many of the features 
observed in these systems. For instance, the shape of 



the KOf-54 brightenings occur for systems at low orbital 
inclinations while larger inclinations can cause the light 
curve to first increase in brightness and then decrease 
(KIC 3547874), or vice versa (KIC 9790355), depending 
on the angle of periastron. Eccentricity has the effect of 
changing the duration of the heartbeat event; heartbeats 
in more eccentric orbits span a shorter phase interval. 

Tidal distortion is not the only physical effect that 
can change the observed brightness of an eccentric bi- 
nary system. Irra d iation, tidally-induced pulsations 
(IKumar et all [l995l: iWillems fc AeTti [200l. Do ppler 
boosting (iBloemen et al.ll2011l : lvan Kerkwiik et al.|[2010l: 
iLoeb fc Gaudill2003f) . and eclipses must also be consid- 
ered in order to obtain the best fit to the light curves. 
Fitting an extensive set of physical effects to each systern, 
as do ne for KOI-54 (Welsh et al. 2011; Burkart et al] 
l2012f ). is complex and beyond the scope of this discovery 
paper. We do, however, show that the tidal distortions 
alone can fit the dominant features of most of the light 
curves. By doing so we also approximate the orbital pa- 
rameters of these systems. 

3.1. Fitting Eccentric Tidal Distortions 

We use the analytic model of lKumar et al.l (|1995l . equa- 
tion 44) for the tidal oscillations induced in eccentric bi- 
naries to obtain the orbital parameters for each of the 
heartbeat stars. Our fit to the fractional flux (^) con- 
tains six parameters: eccentricity (e), orbital inclination 
(i), angle of periastron (w), an amplitude scaling factor 
(S), a fractional flux offset (C), and time of periastron. 
The period and the semi-major axis (a) were fixed to one 
when performing the fit to the folded light curves. R{t) 
describes the distance between the two stars as a func- 
tion of time and (j){t) gives the phase of the ellipse as a 
function of time. 



SF 1-3 sin^ (i) sin^ {4>{t) - uj) 

"f" {R{t)/ar 



C 



(1) 



Inclination is defined such that the plane of the orbit 
is perpendicular to the line of sight at an inclination of 
zero. Because the stellar distortion is symmetric along 
the axis of the elongation, the angle of periastron is re- 
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Figure 3. Folded light curves plotted as relative flux against orbital phase; the zero phase relates to the brightest or dimmest event in 
the folded light curve. Phase in fractional units is labeled along the bottom axis and time in days is labeled along the top axis. The KIC 
number and period are labeled on each plot. We smooth each folded light curve by performing a 15-min running average. 
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0.5 1 1.5 2 

Figure 4. Same as Figurc|3]for the remaining iieartbeat stars. 

ported from —90° to +90°; uj = when looking down 
the minor axis of the ellipse and increases in the direc- 
tion of the orbiting primary. The fit was performed using 
a Levenberg-Marquardt chi-square minimization. 

The results of our fits can be found in Figures [S] and 
[7] and Table [1] Since we fit folded light curves with ex- 
tremely small uncertainties, and did not fit all the fea- 
tures present in the light curve, our fits result in large 
reduced chi-square values with small error bars. To com- 
pensate, the Icr uncertainties on the fitted parameters 
presented in Table [T] represent the error on the fitted pa- 
rameters when the uncertainties on the folded light curve 
are increased so that the reduced chi-square is equal to 
one. The error bars vary greatly depending on the chi- 
square of the fit. 

We have not included several physical parameters 



known to affect the light curves of such highly eccen- 
tric systems. These effects will have a different impact 
depending on the properties of the stars in the system. 
In many cases the harmonic pulsations are large and may 
significantly alter the results of the fit. As an example, 
KIC 4248941 shows large harmonic pulsations that could 
easily sway the results of our fits. Also, the exclusion of 
orbital effects such as irradiation (and to a lesser de- 
gree Doppler boosting) can prevent us from obtaining a 
good fit for the true orbital parameters. For instance 
KIC 4372379, KIC 4936180, and KIC 7622059, all show 
a heartbeat shape but are not fit well by our simple tidal 
distortion model. Adding the additional physical effects 
may account for the discrepancy. 

Wh ile incomplete, we note that the IKumar et al.l 
(|1995f) model can give a good first estimate of the oribtal 



Thompson et al. 



i=70°co=75° 



e=0.4 

i=70°a)=50° 



i=70°(o=25° 



i=70°co=-25'^ 



e=0.6 

i=70°co=-50° i=70°(o=-75° 




0.8 1 1.2 
i=55°a)=75° 



0.8 1 1.2 
i=55°co=50° 



0.8 1 1.2 
i=55°(0=25° 



0.8 1 1.2 
i=55°(0=-25° 



0.8 1 1.2 
i=55°(o=-50° 



0.8 1 1.2 
i=55°(0=-75° 



CD 
T3 

H — ' 

"q. 

E 
CO 



o 




1 


A 


0.5 


/ 





J\r 


-0.5 


V 



1 


A 


0.5 







-\^ 


-0.5 


V 




0.8 1 1.2 
i=40°(o=75° 



0.8 1 1.2 
i=40°co=50° 



0.8 1 1.2 
i=40°co=25° 



0.8 1 1.2 
i=40°co=-25° 



0.8 1 1.2 
i=40°(o=-50° 



0.8 1 1.2 
i=40°(o=-75° 




1 


A 


0.5 


^V 





^ 


-0.5 


V 



0.8 1 1.2 
i=25°co=75° 



0.8 1 1.2 
i=25°co=50° 



0.8 1 1.2 
i=25°(0=25° 



0.8 1 1.2 
i=25°co=-25° 



0.8 1 1.2 
i=25°co=-50° 



0.8 1 1.2 
i=25°(0=-75° 




0.8 1 1.2 0.8 

orbital phase 



Figure 5. Synthetic light curves generated from the analytic model of IKumar et al.l l|1995l ) for tidal distortions in an eccentric binary. 
The inclination (i) increases vertically and periastron angle (oj) increase horizontally. The relative flux of each curve has been normalized 
by the largest change in brightness. Since the grid is symmetric around a; = 0, we present the model at two eccentricities, 0.6 and 0.4 (red 
and blue lines respectively). 



parameters even in cases where the irradiation is sig- 
nificant, and has the advantage of not requiring knowl- 
edge about the properties of the stars to make this es- 
timate. For KOI-54, where the irradiation makes-up 
almost half of the change in flux, the fit to the tidal 
distortion model gives a reasonable estimate of the or- 
bital parameters. According to the detailed modeling of 
iWelsh et al.l (|2011[ ). the orbital parameters of KOI-54 are 
e =0.8315(32), w =39°.46(51), and i =5°. 52(10). Our fit 
is presented in the last panel of Figure [T] While our 
model does not include irradiation, the pulsations, nor 
the information from radial velocity measurements, our 
eccentricity, 0.82(01), angle of periastron, 52°(30), and 
inclination, 11°(3), all match within the two sigma error 
bars. 

The ability to fit KOI-54 well with only the tidal dis- 
tortion model is in part because of its inclination. At 
low inclinations the shape of the contribution from irra- 
diation closely resembles that of the tidal distortion. For 
higher inclinations this will not always be the case, and 



so when the irradiation is significa nt it can appreciabl y 
alter the shape of the heartbeat (see lBurkart et aLll2012f) . 
Once knowledge of the companion can be obtained, the 
effects of irradiation can be included accurately in the 
model of the light curve. 

We have three obviously eclipsing systems: 
KIC 3230227, KIC 4936180, and KIC 8719324. We fit 
these systems with the inclinations of 42.79°, 23.2°, 
and 75.78°, respectively. Our model for KIC 4936180 
is obviously incomplete and so the fitted inclination is 
likely incorrect. Both KIC 8719324 and KIC 3230227 
have very good fits, but rather small inclinations for 
eclipsing systems. It is possible to achieve an eclipse 
at these inclinations if the distance at periastron for 
these highly eccentric systems is very small (^ SR©). 
This is possible if the secondary has a rather low mass 
(^ 0.5 Mq). a more likely possibility is that these are 
cases where the irradiation has a significant effect on 
the shape of the light curve, causing us to estimate the 
inclination at a lower value. 
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This simple model, despite being incomplete, allows 
us to illustrate that tidal distortions in an eccentric bi- 
nary system can explain many of the features in the 
light curves. Fitting these features will allow us to es- 
timate the orbital parameters of these systems. One 
trend noticed from fitting our simple model is that the 
eccentricity is smaller for shorter period systems; see 
Figure ID While we are less likely to find long-period, 
low-eccentricity systems with observable tidally-induced 
variations, we would expect to detect highly eccentric, 
short-period systems if they existed. Since shorter pe- 
riod orbits undergo stronger tidal forces, they are likely 
to circularize more quickly, leaving a dearth of highly ec- 
centric short period bina ries. Similar trends have been 
seen bv lPrsaet all ()2008[ ) and lMathieu fc Mazehl (|1988[ ) 
for eclipsing binaries. 

4. SPECTRA 

In order to measure the stellar parameters and to 
confirm the prediction that these are eccentric binary 
systems, we obtained optical spectra of the 17 differ- 
ent heartbeat stars at the KPNO Mayall 4-m telescope 
between 2011 May 14 and September 7 (UTC) using 
the long-slit spectrograph RCSpec and T2KA Camera. 
The KPB-22b grating was used in 2nd order with a slit 
width of 1.0" to produce a dispersion of 0.72Apixel~^ 
and resolution R ~ 4000. Each spectrum was observed 
with an exposure time between 180 and 900 s, depending 
on brightness and conditions and each object exposure 
was preceded by a 30-s exposure of an FeAr arc lamp 
spectrum for wavelength calibration. A standard star 
(one of BD-H40 4032, Feige 34, or Wolf 13 46) was ob- 
served each night usi ng the same method (jStond ITQTTI : 
iBarnes &: Havesiri982[ ). Bias exposures and quartz lamp 
flat fields were taken the day before each night of obser- 
vations for calibration. 

We reduced the spectra using standard IRAF (Image 
Reduction and Analysis Facility) tools. First, the over- 
scan bias level was subtracted from each frame. Next, 
a master residual bias and normalized flat-field frame 
were constructed from sets of bias and quartz lamp ex- 
posures and then applied to each spectrum frame to cor- 
rect for any remaining 2-D bias patterns and pixel-to- 
pixel variations in detector response. The IRAF oned- 
spec package was used to define apertures and to extract 
sky-subtracted spectra. The same apertures are also em- 
ployed to extract an arc lamp spectrum for each object. 
The spectra are then resampled on a linear wavelength 
scale using fits to the arc lamp spectra. The reduced 
standards were then fit to the flux standards provided 
with the IRAF onedspec package to produce a sensitiv- 
ity curve in order to correct each science spectrum. 

These spectra suffer from degraded focus outside of 
the range 3800 — 4700 A where the fluxes are difficult to 
properly calibrate. Consequently, we have restricted our 
spectra fits to this limited wavelength range. 

We have modeled each spectrum usin g a library of 
synthe tic model atmosphere spectra from iMunari et all 
(I2005D. which are based on the model atmospheres of 
iCastelli fc Kuruc3 ()2003l ). The predicted spectra are re- 
sampled to a resolution to match the observed spectra 
and the best-fitting model, scaled in flux to best match 
the spectrum, is found by minimizing the rms residuals 
in the fit over the wavelength interval 3800 — 4700 A. The 



grid of atmosphere models are spaced at discrete inter- 
vals on Tcff (250K) and log5(0.5 cgs), and so the best 
fits have values at these grid points with uncertainties at 
least as large as half of the grid point spacing. See Fig- 
ure [9] for individual spectra of each system over-plotted 
with the best fit model. The results are listed in Table [TJ 

Our spectral fits assume that each target is dominated 
by the primary in the system. In all but one case, we only 
see spectral lines for one star, even when taking obser- 
vations near periastron. The exception is KIC 3547874 
which shows evidence of being a double lined binary like 
KOI-54 (see Figure [9]) . The spectrum that reveals the 
double lines was taken at a phase concurrent with the 
heartbeat. From the separation in the spectral lines, the 
spectrum indicates that the primary has a radial veloc- 
ity of lOOill kms~^ relative to a spectrum taken just 
before the heartbeat. 

Since our spectra are of moderate resolution, several 
systems may show evidence of double lines if spectra 
were taken at the correct phase with a higher resolu- 
tion spectrograph. However, in order for the secondary 
to contribute a significant fiux, it would have to have a 
similar temperature as the primary (or be a giant star). 
A blend of two stars with similar temperatures will not 
cause a large error in our fitted temperatures and gravi- 
ties. Note that KIC 8719324 and KIC 4248941 both had 
spectra taken when the radial velocity shift was much 
larger than the spectral resolution and therefore cannot 
be double line binaries. Both would have shown evidence 
of double lines if the companion was of a similar temper- 
ature. 

4.1. Radial Velocity Measurements 

For five systems, we obtained more than one spectrum 
and were able to measure the relative velocities at differ- 
ent orbital phases. We measured the radial velocity be- 
tween the spectral lines by performing a cross-correlation 
against one spectrum. We show relative radial velocity 
measurements for these systems phased with the folded 
light curve in Figure [101 AH radial velocities have been 
corrected for the Earth's motion around the sun. The er- 
ror bars are estimated from the full width half maximum 
of a Gaussian fit to the peak of the cross-correlation func- 
tion. Note that KIC 3547874 is a double- lined binary and 
so our cross-correlation technique is too simple to accu- 
rately measure the radial velocities, especially when the 
lines are not well resolved. 

4.2. Fitting Radial Velocity Measurements 

For four systems we have a time series of radial velocity 
measurements that show a distinct non-sinusoidal shape, 
indicative of an eccentric binary system; see Figure 1101 
Unfortunately we do not have enough data for any of 
these systems to constrain the orbital parameters of the 
system; however we can show that the radial velocity 
measurements are consistent with the orbital parameters 
determined from the fits to the folded light curve. We 
fit the radial velocity curve constraining the eccentricity, 
angle of periastron, and time of periastron passage to 
those listed in Table [1] The velocity semi- amplitude was 
a free parameter of the fit. See the blue line in Figure [TU] 
for the results of this fit. The parameters of the fit are 
labeled on the figure for each system. 
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Figure 6. Best fits of the analytic model of lKumar et alj l|1995l ) (dashed green line) to folded light curves (solid black). The eccentricity, 
inclination, and angle of periastron (e, i, and co respectively) established from the fit are noted at the bottom of each plot. The lower panel 
for each set contains the residuals of the fit. Our model excludes many known second-order physical effects, leaving large variations in the 
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For KIC 3230227, KIC 4248941, KIC 4372379 and 
KIC 8719324 the tidal distortion model is mostly consis- 
tent with the radial velocity fits. Those remaining incon- 
sistencies are in part caused by using orbital parameters 
determined by fitting an incomplete model to the light 
curves. However, another possible reason for the incon- 
sistency is that the radial velocity measurements may 
include the unmodelled tidally-induce d distortions and 
oscillations on the surface of the star. iWillems fc Aerta 
([2002) discuss these velocity fields and show that for high 
inclinations, the amplitude of these velocities may be as 
high as 15 km s'"^ for one pulsation mode on a B-type star 
with an eccentricity of 0.5. Despite being cooler, signifi- 



cant pulsation velocity fields may exist on our heartbeat 
stars. Factors that increase the size of the tidal force, 
such as eccentricity and the mass of the secondary, can 
also increase the size of this effect. Note, this would 
not have been a dominant effect for KOI-54 because of 
the low inclination of that system. Whether we are see- 
ing these surface velocity fields is unclear without bet- 
ter time-resolved radial velocity measurements, but such 
considerations must be made when fitting the data on 
these complicated, tidally dynamic systems. 



5. DISCUSSION 
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metric tidal distoriton model plotted against the orbital period for 
the 17 Kepler s ystems. KOI-54 is plotted as a 'x' and HD 174884 
liVlaceroni et al.| ([2009) is plotted as a square. The error bars are 
smaller than the size of the points. 



Our collection of heartbeat stars were discovered 
serendipitously while searching for planets in the Ke- 
pler archive. Kepler targets were intentionally selected to 
look for "Sun-like" stars, and so primarily G and K-type 
stars were selected. For stars brighter than 15.5 mag, 
97 per cent of the selected stars are coole r than 6500 K 
with the distribution peaking near 5500 K (jBatalha et al.l 
120101 ). The temperatures of the observed stars in our 
group of binaries range from 6000 to 9250 K. While we 
certainly have not found all eccentric systems with dy- 
namic tidal distortions in the exoplanet data, we might 
expect our collection to be dominated by G-type stars if 
heartbeat stars were common across all main sequence 
stars. 

The dominance of hot stars in our sample is likely the 
result of the size of the tides on the primary star. The 
fractional amplitude of the tidal distortions depends on 
the star's surface gravity a s well as the size of the tidal 
force from the secondary (jPfahl et al.ll2008D . An A5V 
star has a lower log g and will have approximately twice 



the amplitude of a G5V star undergoing the same tidal 
force, making it easier to detect. We may be seeing this 
trend in Figure [TT] where we plot amplitude against ef- 
fective temperature. We measured the amplitude as the 
absolute value of the largest deviation from zero level in 
the folded light curves (without including the amplitude 
of the obvious eclipsing events). We have not adjusted 
this plot for effects that can affect the strength of the ex- 
ternal tidal forces (such as eccentricity and the mass of 
the companion relative to the primary). We do note that 
no obvious trend between temperature and eccentricity 
exists in this collection of heartbeat stars (see Table [IJ . 
As we continue to search the Kepler data set and can 
fold longer data sets, we will likely find a larger popu- 
lation of cooler main sequence heartbeat stars with low 
amplitudes. 

5.1. Circularization 

From the perspective of stellar evolution, highly eccen- 
tric short-period binary systems are a mystery. Highly 
eccentric systems with small periastron separations are 
expected to tend towards circular orbits at a rapid rate 
as their tidal interactions are efficient at dissipating en- 
ergy (tidal interactions scale as (ii!*/a)~^). Thus the lack 
of short-period, highly-eccentric stars in Figure [8] is not 
unexpected. 

There are currently two theories that describe the cir- 
cularization of binary orbits: the equilibriur n tide mode l 
(lHutlll980D and the d ynamical tide model (|Zahnlll975D . 
lEggleton et al.l ()1998[ ) states that the former is most 
likely applicable to stars that have relatively small radii 
with respect to their periastron distance and the latter 
for systems where the stellar surfaces are in close proxim- 
ity at periastron. The equilibrium tide model attributes 
orbital circularization to the torque caused by the mis- 
alignment of the tidal bulges on the stellar surface with 
respect to the star's instantaneous equipotential shape. 
This model suggests that in certain cases, when the ratio 
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Table 1 

Eccentric Binary Systems 



KIC 


2MAj^9 


^A.ag 


quarters 


■Ks 


log 9 


t-eriod (d) 


e 


u (") 


^(") 


3230227 


19202701+3823595 


9.002 


0-3 


8750 


5.0 


7.04711(87) 


0.5880(40) 


-67.9(1.2) 


42.79(46) 


3547874 


19290789+3839565 


10.952 


0-6 


6750 


4.5 


19.69284(59) 


0.6480(20) 


40.40(20) 


46.434(56) 


3749404 


19281908+3850135 


10.575 


0-6 


7500 


5.0 


20.306630(60) 


0.6360(20) 


-58.17(75) 


37.53(17) 


4248941 


19102476+3923516 


12.159 


0,2-7 


6750 


4.5 


8.644718(34) 


0.423(14) 


-50.5(4.3) 


68.3(5.5) 


4372379 


19340487+3926520 


13.810 


1-5,7 


6750 


4.5 


4.535167(35) 


0.342(92) 


77(26) 


48(11) 


4847343 


19394839+3955460 


12.837 


0-6 


6250 


4.5 


11.4162193(59) 


0.4050(30) 


-42.81(62) 


24.23(16) 


4936180 


19402728+4004162 


15.270 


1-6 


6750 


5.0 


4.640867(59) 


0.079(58) 


8(16) 


23.2(6.7) 


5034333 


19493609+4008581 


11.520 


1-7 


9250 


4.5 


15.29434(21) 


0.5750(20) 


-81.41(38) 


48.21(13) 


7622059 


19434079+4316543 


13.272 


1-7 


6250 


4.5 


10.403145(43) 


0.406(91) 


-63(16) 


28.0(5.5) 


8264510 


20035245+4411313 


10.609 


0-7 


8250 


4.0 


5.6867311(64) 


0.3380(30) 


-60.42(49) 


29.03(16) 


8719324 


20042742+4449328 


11.614 


0-6 


7750 


4.5 


10.23262347(27) 


0.6000(20) 


-16.61(42) 


75.78(89) 


9790355 


19545806+4634382 


14.055 


3-6 


6000 


3.0 


14.56034(73) 


0.5130(70) 


-55.1(2.3) 


46.33(69) 


9835416 


19344716+4637138 


13.338 


4-6 


7250 


5.0 


4.03650(11) 


0.279(23) 


-44.3(7.5) 


49.8(3.5) 


9899216 


19403879+4645023 


10.886 


1-6 


7500 


4.0 


10.915763(20) 


0.6470(40) 


81.2(1.3) 


20.74(61) 


10873904 


19462402+4817456 


13.065 


1-6 


6250 


5.0 


9.88553(23) 


0.4360(20) 


-12.67(18) 


42.196(47) 


11494130 


18560643+4924555 


10.988 


0-3 


6750 


4.5 


18.9555(14) 


0.6180(40) 


62.3(1.6) 


48.29(52) 


11568657 


19450715+4931040 


13.246 


0-4 


6000 


4.5 


13.47697(98) 


0.5650(20) 


-46.80(30) 


48.202(96) 



Note - The values for T^ff and logg are the result of fits to a grid of model spectra with a spacing of 250 K and 0.5 (cgs), respectively 
(seejffl. The eccentricity (e), angle of periastron (uj), and inclination (i) are the result of the fits to the tidal distortion model described 
in i|3.1l This model is likely incomplete for many of these systems; use the resulting orbital parameters with caution. The error bars 
represent la uncertainties and do not include the error associated with any systematic effects. 
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of the spin to the orbital angular momentum is large, the 
orbital evolution of the system can be erratic; eventually 
leading to circularization, collision or even escape. 

For stars that interact more dramatically during peri- 
astron (as li kely th e case for these highly eccentric heart- 
beat stars). IZahnl (J1975 ) has shown that the dynamical 
tidal forces generate oscillation in the surface layers of the 
star that are damped through tidal interactio ns; a trend 
confirmed bv lKhaliulhn fc KhahuUinal (|2011f ). This pro- 
cess is highly effective at energy loss and results in a 
circularization time scale much shorter than the lifetime 
of ho t main sequence stars (on the order of 10^ y, iHutI 
[1981 . 

One way to maintain eccentric binaries in hot main se- 
quence stars is by including a third bo dy. Advances on 
the w ork of iMavor fc Mazehl ()1987[ ) by iTokovinin et al.l 
(|2006!) recently determined that in binary systems with 
periods less than 3 d the frequency of tertiary compo- 
nents is as large as 96% (dropping to 34% for binary 
systems with periods > 12d). A tertiary component can 
both remove angular momentum from and perturb the 
inner binary system, allowing it to maintain a high ec- 
centricity despite the strong dynamic tidal forces. Thus 
the existence of highly eccentric binary systems, such as 
those in our sample, is possibly a consequence of the 
presence of a third body. 

5.2. Induced Pulsations 

The tidal forces also induce pulsations on these stars. 
We do no analysis of these obvious features in the light 
curve, but we do note their presence. A cursory analy- 
sis of the Fourier transform and folded light curves show 
that a majority have pulsations that are harmonics of 
the orbital period and several also show modes indepen- 
dent of the orbital period. The pulsations with a period 
at an integer harmonic of the orbit al period are likely 
driven by the variable tidal forc es (jWelsh et al.l 120111 : 
lWillemsll200l iKumar et aP 119951 ). Other modes could 
be caused by one star being a normally pulsating star, 
or, the tidal forces are exciting additional m odes through 
nonlinear couplings with re sonant modes (jFuUer fc Lail 
IMllWufc Goldreich|[200l . 

Since the tidally-induced modes are harmonics of the 
orbital frequency, the only way to accurately measure the 
pulsati ons is to first corre ctly fit and remove the orbital 
effects (Welsh et al.ll201ll ). Because of the careful mod- 
eling effort required for each individual star, we leave the 
extraction and analysis of the tidally-induced pulsations 
for subsequent papers. 

6. CONCLUSIONS 

With the new level of precision and coverage afforded 
by the Kepler time series data, we have discovered 17 
eccentric, tidally distorted, binary stars we have called 
heartbeat stars. KOI-54 is a long period, highly eccen- 
tric, example of this class. The variety of light curve 
shapes arises primarily from variable tidal distortions be- 
ing observed from different inclinations and periastron 
angles. Radial velocity measurements, as well as fits to 
the light curves, support this model. As Kepler obtains 
more data, we will likely find more examples of stars in 
this class with lower amplitudes. 

Detailed modeling of these systems will be challenging. 
The tidal forces are changing throughout the orbit caus- 



ing large tidal distortions and induced pulsations on pos- 
sibly both stars. These effects can cause non-sinusoidal 
variations in both the light curve and the radial velocity 
measurements. Plus, large orbital velocities can cause 
Doppler boosting to add to the ffux variations while 
the secondary can cause significant irradiation variations. 
When more complete radial velocity measurements are 
obtained, the parameters of the companion can be con- 
strained and a better model can be developed for each 
star. With detailed models, these stars offer unique lab- 
oratories in which higher order effects can be studied, 
effects often easily ignored in stellar structure calcula- 
tions. 

This is the first opportunity to observationally study 
dynamic tidal interactions on a class of stars. Our first 
analysis of each system has already shown some of the 
potential of these systems to constrain theories regard- 
ing energy transfer in tidally dynamic systems. From our 
fits to these heartbeat systems, we confirm that highly 
eccentric binaries with very short periods must be rare, 
indicating that tidal friction may have circularized these 
systems. Further investigation may confirm the theories 
of Zahn (1975) and KhahuUin & Khaliulhna (2011) re- 
garding the time scales of circularization. Additionally, 
the strength and period of the tidally-induced pulsations 
varies across these systems. With a class of such systems, 
we can begin to decipher how this driving mechanism 
excites certain oscillation modes in a diverse sample of 
stars. 
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Figure 9. Spectrum of each object in black. Model fits to the spectra are over-plotted in red. Physical parameters of the model are 
specified in the upper right hand corner of the plot. The last panel shows the double-lined spectrum of KIC 3547874 obtained at an orbital 
phase of 0.986. 



Eccentric Binaries with Dynamic Tidal Distortions 



13 




KIC 8719324 
P=10.2326 d 



I I 




"A=48.9 km/s e = 0.60 w = -16.6 



1.5 



----<)- - 



KIC 3230227 
p= 7.0471 d 




7.8 km/s e = 0.59 w = -67.9 



1 1.5 
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eccentricity, angle of periastron and time of periastron passage to the values measured in i|4.2l (blue line). The amplitude {A), eccentricity 
(e) and angle of periastron (w) of the radial velocity fit are listed on the plot. 
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Figure 11. The observed amplitude of the heartbeat plotted 
against the measured temperature of the primary star. The error 
in the amplitude is approximately the size of the point. 
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